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ABSTRACT 

Context. In April 2013, the nearby (z=0.031) TeV blazar, Mkn 421, showed one of the largest flares in X-rays since the past decade. 
Aims. To study all multiwavelength data available during MID 56392 to 56403, with special emphasis on X-ray data, and understand 
the underlying particle energy distribution. 

Methods. We study the correlations between the UV and gamma bands with the X-ray band using the z-transformed discrete correla¬ 
tion function. We model the underlying particle energy spectrum with a single population of electrons emitting synchrotron radiation, 
and do a statistical fitting of the simultaneous, time-resolved data from the 5w!/t-XRT and the AmSTAR. 

Results. There was rapid flux variability in the X-ray band, with a minimum doubling timescale of 1.69 ±0.13 hrs. There were no 
corresponding flares in UV and gamma bands. The variability in UV and gamma rays are relatively modest with ~ 8% and 16% 
respectively, and no significant correlation was found with the X-ray light curve. The observed X-ray spectrum shows clear curvature 
which can be fit by a log parabolic spectral form. This is best explained to originate from a log parabolic electron spectrum. However, a 
broken power law or a power law with an exponentially falling electron distribution cannot be ruled out either. Moreover, the excellent 
broadband spectrum from 0.3 - 79 keV allows us to make predictions of the UV flux. We find that this prediction is compatible with 
the observed flux during the low state in X-rays. However, during the X-ray flares, depending on the adopted model, the predicted flux 
is a factor of 2 - 50 smaller than the observed one. This suggests that the X-ray flares are plausibly caused by a separate population 
which does not contribute significantly to the radiation at lower energies. Alternatively, the underlying particle spectrum can be much 
more complex than the ones explored in this work. 
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1. Introduction 

According to the unification s cheme of Active Galactic Nuclei 
(AGNs) by lUrrv & Padovanil (Il995h . blazars are a subclass of 
AGNs with a relativistic jet aligned close to the line of sight. 
Blazars are further subdivided into BL Lacs and FSRQs where 
BL Lacs are characterized by the absence of (or very weak) 
emission lines. They show high optical polarization, intense 
and highly variable non-thermal radiation throughout the entire 
electromagnetic spectra in time scales extending from minutes 
to years, apparent super-luminal motion in radio maps, large 
Doppler factors and beaming effects. The broadband spectral en¬ 
ergy distribution (SED) of blazars is characterized by two peaks, 
one in the IR - X-ray regime, and the second one in y-ray regime. 
According to the location of the first peak, BL Lacs are further 
classified into Low energy peaked BL Lacs (LBLs) and Hig h 
energy peaked BL Lacs (HBLs) (IPadovani & Giommil Il99-5h . 
Both leptonic and hadronic models have been used to explain 
the broadband SED with varying degrees of success. The origin 
of the low energy component is well established to be caused 
by synchrotron emission from relativistic electrons gyrating in 
the magnetic field of the jet. However the physical mechanisms 
responsible for the high energy emission are still under debate. 
It can be produced either via inverse Compton (IC) scattering 
of low frequency photons by the same electrons responsible for 
the synchrotron emission (leptonic models), or via hadronic pro¬ 


cesses initiated by relativistic protons, neutral and charged pion 
decays or muon cascades (hadronic models). The seed photons 
for IC in leptonic models can be either the synchrotron pho¬ 
tons itself (Synchrotron Self Compton, SSC) or from external 
sources such as the Broad Line Region (BLR), the accretion disc, 
the cosmic microwave background, etc (External Comp ton, EC). 
Fora comprehensive review of these mechanisms, see iBottcherl 
(I2007h . 


Mkn 421 is the closest (z = 0.031) and the most well stud¬ 
ied Te V blazar. It vvas als o the first detected extragalactic TeV 
source (IPunch et al.lll992h . and one of the brightest BL Lac ob¬ 
jects seen in the UV and X-ray bands. It is a HBL, with the 
synchrotron spectrum peaking in the X-ray regime. Moreover, 
the X-ray emissi on is known to be highly correlated with 
the TeV emission (Blazeiowski et al.l2005HK. Katarzvnski et all 


l2005t Hie Qian et ^ 1998 ), but like other HBLs, shows mod¬ 


erate correlation with the GeV emission (iLi etal.ll2013h . It i 


IS 


episodes bv several authors (Aleksic et al. 

20121; Shuklaetal. 

2012 

; llsobe et al.l 201(1 Rrawczvnski et al. 

20011 Aleksic et al. 

2010 

; lAcciari et al.l 20091 lUshio et al. 20091 iTramacere et al. 

2009 

; iHoran et al. 

20091 Lichti et al.l 20081; Fossati et al. 2008; 

Albert et al. l2007l 

Brinkmannet al.ll2003). A detailed study of 


(1201 ih with the most well sampled SED till date. 
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During April 2013, Mkn 421 underwe nt one of the large st X- 
ray flares ever recorded in the past decade dPian et al.l2014t) . The 
source was simultaneously observed by Swift and AmSTAR dur¬ 
ing this flaring episode, and we use these observations to study 
the spectral variations. To strengthen our study further, we sup¬ 
plement the X-ray information with other multiwavelength ob¬ 
servations available. 

The main focus of our work being the joint spectral fitting 
between the Nuclear Spectroscopic Telescope Array AmSTAR 
the Swift-XRT telescopes, and using this, we investigate the un¬ 
derlying particle energy spectrum. The high photon statistics 
during the flare, coupled with the excellent spectral response of 
Swift-XRT and AmSTAR gives us a rich X-ray spectrum from 
0.3 -19keV to explore. A study of the qu iescent state of Mkn 42 1 
using AmSTAR data was performed by iBalokovic et all (l2013h . 
In Section |2] we describe the data reduction techniques from the 
various instruments. Section[3]lists the multiwavelength tempo¬ 
ral results, while the X-ray spectral modelling is described in 
Section |4] We discuss the implications of the results in Section 
0 

2. Multiwavelength Observations and Data Analysis 

The huge X-ray flare of Mkn 421 during April 2013(2013 April 
10 to 21; MJD 56392-56403) is simultaneously observed by the 
AmSTAR and Swift X-ray and UV telescopes. The y-ray behav¬ 
ior of this source during this flare was obtained by analyzing 
Fermi-TP^ observations. In addition, we also include the X-ray 
observations by MAXI and optical observations by SPOT for the 
present study. The analysis procedures of these observations are 
described below. 

2.1. Fenni-Large Area Telescope Observations 

The Fenni-LAT data used in this work were collected cov¬ 
ering the period of the X-ray outburst (MJD 56390-56403). 
The standard data analysis procedure as mentioned in the 
Fermi-LAT documentatiorQ was employed. Events belonging 
to the energy range 0.2-300 GeV and SOURCE class were 
used. To select good time intervals, a filter “DATA_QUAL>0” 
&& “LAT_C0NFIG==1” was chosen and only events with less 
than 105° zenith angle were selected to avoid contamination 
from the Earth limb y-rays. The galactic diffuse emission 
component gll_iem_v05_revl.fits and an isotropic component 
iso_source_v05_revl.txt were used as the background models. 
The unbinned likelihood method included in the pylikelihood 
library of Science Tools (v9r33p®) and the post-launch in¬ 
strument response functions P7REP_SOURCE_V15 were used 
for the analysis. All the sources lying within 10° region of inter¬ 
est (ROI) centered at the po sition of Mkn 421 and defined in the 
second Fenni-LAT catalog (iNolan et al.ll2012h . were included in 
the xml file. All the parameters except the scaling factor of the 
sources within the ROI are allowed to vary during the likelihood 
fitting. Eor sources between 10° to 20° from the centre, all pa¬ 
rameters were kept frozen at the default values. The source was 
modelled by a power law as in the 2EGL catalog. 

2.2. NuSTAR Observations 

AmSTAR (iHarrison et al.l 1201^ features the first focussing X- 
ray telescope to extend high sensitivity beyond 10 keV. There 
were 11 AmSTAR pointings between the aforementioned dates, 

' http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 


the details of which are given in Table [T] The AmSTAR data 
were processed with the NuSTARDAS software package v. 1.4.1 
available within HEASOET package (6.16). The latest CALDB 
(v.20140414) was used. After running nupipeline v.Q.4.3 
on each observation, nuproducts v.0.2.8 was used to ob¬ 
tain the lightcurves and spectra. Circular regions of 12 pix¬ 
els centered on Mkn 421, and of 40 pixels centered at 165.96, 
38.17 were used as source and background regions respec¬ 
tively. The spectra from the two detectors A and B were 
combined using addascaspec, and then grouped (using the 
tool grppha V.3.Q.1) to ensure a minimum of 30 counts 
in each bin. To get strict simultaneity with Swift-XRT ob¬ 
servations, observation id 60002023025 was broken into 4 
parts 56393.15591890 to 56393.29538714, 56393.29538714 
to 56393.91093765, 56393.91093765 to 56393.96788711 and 
56393.96788711 to 56394.37822500 MJD. 

2.3. Swift Observations 

There were 15 Swift pointings between the aforementioned 
dates, the details of which are given in Table |2l Publicly avail¬ 
able daily binned source counts were taken from the 5w//f-BAT 

webpagq3- 

The XRT data (iBurrows et al.l 120051) were processed with 
the XRTDAS software package (v.3.0.0) available within 
HEASOET package (6.16). Event files were cleaned and cal¬ 
ibrated using standard procedures (xrtpipeline v.Q. 13.®), 
and xrtproducts v. ®. 4.2 was used to obtain the lightcurves 
and spectra. Standard grade selections of 0-12 in the Windowed 
Timing (WT) mode are used. Circular regions of 20 pixels cen¬ 
tered on Mkn 421 (at 166.113 and Dec 38.208) and of 40 pix¬ 
els centered at 166.15, 38.17 were used as source and back¬ 
ground regions respective ly. Eor the observati ons affected by 
pileup (counts > lOOc/s) (iRomano et aLlBOOg) . an annular re¬ 
gion of inner radius 2 pixels and outer 20 pixels was taken as 
the source region. The lightcurves were finally corrected for 
telescope vignetting and PSE losses with the tool xrtlccorr 
V . ®. 3.8. The spectra were grouped to ensure a minimum of 30 
counts in each bin by using the tool grppha v. 3 . ®. 1. 

5w!/f-UVOT (iRoming et al.l200^ operated in imaging mode 
during this period, and for most of the observations, cy¬ 
cled through the UV filters UWl, UW2 and UM2. The tool 
uvotsource v.3.3 was used to extract the fluxes from each 
of the images using aperture photometry. The observed magni¬ 
tudes were corrected fo r galactic extinc t ion (E b-v - 0.019 mag) 
using the dust maps of ISchlegel et al.l (Il998l) and converted to 
fl ux units using the zero point magnitudes and conversion factors 
of iBreeveld et al.l (1201 ll) . The tool £lx2xsp v. 2.1 was used to 
convert the fluxes to pha files for use in XSPEC. 

2.4. Other Multiwavelength data 

Publicly available daily binned source counts were plotted for 
MAXI Q As a part of the Fermi multiwavelength support pro¬ 
gram, the SPOT CCD Imaging/Spectro polarimeter at Ste ward 
Observatory at the University of Arizona temith et al.l2009h reg¬ 
ularly observes Mkn 421. The publicly available optical V-band 
photometric and linear polarization data were downloaded from 
their website Q. 


^ http://swift.gsfc.nasa.gov/results/bs70mon/SWIFT_J 1104.4p3812 
^ http://maxi.riken.jp/ 
http://james.as.arizona.edu/~psmith/Fermi/ 
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3. Multiwavelength Temporal Study 


4. X-ray spectral analysis 


The Multi wave length lightcurve during the 10 days period, from 
optical to gamma ray energies, along with the optical polariza¬ 
tion measurements, is plotted in Fig. [1] While there were two 
huge flares in X-rays (on MJD 56395 and 56397), where the 
flux went up by a factor of 10 , the fluxes in the other bands were 
not very variable on the timescale of days. We compute the z- 
transformed discrete correlation using a freely available Fortran 
90 code w i th the details of the method employed described in 
lAlexand^ (Il997l) . We find no lag between the soft (0.3 -10 keV, 
Swift-XRT) and the hard X-ray (3.0 - 79 keV AmSTAR) bands. 
There is no correlation seen between the UV flux and the X-ray 
flux, (zdcfmax - 0.62 + 2.3, at a lag of 2.2 days). Also, while 
the UV flux does not show correlation with the optical polariza¬ 
tion {zdcfmax - 0.61 + 1.1, at a lag of 2.4 days), the X-ray flux 
shows a tighter correlation {zdcfmax - 0.81 + 0.6, at a lag of 3.5 
days) with the latter. There was also a large change in the angle 
of polarization during the two X-ray flares. 

The hardness ratios (computed here as the ratio between 
the 10 - 79 keV co unt rate and the 3-10 keV count rate 
(iTomsick et alJ 1201^ 1 are plotted in Fig. |2] A trend of spec¬ 
tral hardening with increasing flux (Spearman rank correla¬ 
tion, rs=0.58, p - 5.0 ■ 10“*), is o bserved and the same 
is often reported for this source (eg: iBalokovic et al.l 12013t 
IW. Brinkmann et ^l2001h . Moreover, the correlation is much 
tighter during the rising part of the 2 flares (rs=0.92, p - 
4.0 ■ 10“^). These interesting features advocate us to perform a 
more detailed spectral study, which we describe in SectionlH 

The fractional variability amplitu de parameter F^ar 
dVaughan et al.l 1200.31 : IChitnis et al.l l2009l) . computed on daily 
timescales, is used to quantify the multi-wavelength variability. 
It is calculated as 


The close simultaneity between Swift-XRT and VmSTAR ob¬ 
servations allows us to do a joint spectral fitting using XSPEC 
package v. 12. 8 .2. The time periods which have been fitted 
together are shown in Fig. |3] The bin widths are selected as one 
bin per Swift-XRT observation (except for Obs. id. 00035014062 
and 00035014069, which last only for a few minutes) lead¬ 
ing to 13 time bins which we denote as fl - fl3. The state D 
has no Swift-XRT data, whereas the state fl3 has no AmSTAR 
data. Again in Table |3] we list the Swift-XRT and corresponding 
NuSTAR data which have been combined. 

While fitting the broadband X-ray spectrum (0.3 - 79 keV), 
the XRT and the NuSTAR spectral parameters were tied to each 
other, except the relative normalization between the two instru¬ 
ments. To correct for the line of sight absorption of soft X-rays 
due to the interstellar gas, the neutr al hydrogen column density 
was fixed at Nh - 1.92 x lO^^’cm^^ (iKalberla et al.ll200^ . 

4. 1. Fitting the photon spectrum 

It is known that the X-ray spectrum of Mkn 421 shows s ignif- 
icant curvature (iFossati et al.l 120005 iMassaro et al.llT004t) . and 
consistently we also noted that the data cannot be fitted satis¬ 
factorily by a simple power law. On the other hand, a power 
law with an exponential cutoff gives a much steeper curvature 
than observed, yielding to unacceptable fits. A sharp broken 
power law also gives large values in most cases, which sug- 
gests a smooth intrinsic curvature in the spec trum. So, following 
(IMassaro et al.ll2004l iTramacere et alJl2007l) we fit the observed 
spectrum with a log parabola given by 

dN/dE = ( 4 ) 


F 


var ~ 



( 1 ) 


where is the mean square error, x the unweighted sample 
mean, and S'^ the sample variance. The error on F^ar is given as 


O'A.r 


A 


2A x^F,, 



( 2 ) 


Here, N is the number of points. 

The variability amplitude is the maximum for the X-ray 
bands {F^ar ~ 0.75 + 0.10), and significantly lower for the UV 
bands {Fyar ~ 0.08 + 0 . 02 ), suggesting that the emission may 
probably arise from different components in the two bands. The 
variability in the GeV range is also small {F^ar ~ 0.16 + 0.07). 
This goes against the general trend found in blazars that F^ar is 
the maximum for the y-ray band and decreases with frequency 
(IZhang et al.ll2005l : lPaliva et al.ll2015[) . 

We scan the Swift-XRT and the AmSTAR lightcurves for the 
shortest flux doublin g timescale using the following equation 
dFoschini et al.ll 20 Tlh 

F{t) = F{to).2^‘-'«^/^ (3) 


where gives the spectral index at Eb^s- The point of maximum 
curvature, Ep^s is given by 

Ep,, = £4U0(2-“9/2A (5) 

During fitting, £ 4 ^ is fixed at IkeV. In Table |3l we give the re¬ 
sulting reduced 4 f^s for the case of the broken power-law and the 
log parabola models; while in Table |4] we give the fit parameters 
corresponding to the latter model. 

Interestingly we noticed a strong anti-correlation between 
and flux {rs - -0.98, p < 2.2 ■ 10“*®), and a strong correlation 
between flux and Ep^s {rs - 0.97, p < 2.2 • 10“*®) which implies 
that during flares, the spectral index at IkeV hardens and the 
peak of the spectrum shifts to highe r energies. This behaviour o f 
the source has often been reported (IMassaro et al.ll2004l l2008li . 
In addition, we did not see any correlation betw een a, and Es 

i rs =■ -0.35, p - 0.24), which was seen by IMassaro et al.l 
20041) . We also noticed that there is no correlation observed be¬ 
tween the curvature parameter (3 and the peak of the curvature 
Ep {rs = 0.63, p - 0.02). These cross plots are shown in Figure 


4.2. Emitting Particie Distribution 


where F{t) and £(fo) are the fluxes at time t and to respectively, 
and T is the characteristic doubling/halving time scale. The 
fastest observed variability in the AmSTAR band is 1.69 + 0.08 
hours between MJD 97.06438 to 97.08573. This is compara¬ 
ble to the ve ry fast variability ob served in this source with the 
Beppo-SAX (iFossati et al.ll2000^ . 

The above study is performed only for those periods where 
the flux difference is significant, at least at the 3cr level. 


The excellent spectral resolution of AmSTAR gives us an un¬ 
precedented view of the high energy X-ray behavior beyond 
2QkeV. Coupled with Swift-XRT, we have, for the first time, an 
uninterrupted, well resolved spectrum from 0.3 - 79 keV. This 
allows us to go beyond only fitting the photon spectrum with 
various spectral forms. Rather, in this work, we try to study the 
underlying particle distributions which give rise to the observed 
photon spectrum. 
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We consider the case where X-ray emission arise from a rel¬ 
ativistic distribution of electrons emitting synchrotron radiation. 
The electrons are conhned within a spherical zone of radius R 
hlled with a tangled magnetic field B. Due to relativistic motion 
of the jet, the radiation is boosted along our line of sight by a 
Doppler factor 6. A good sampling of the entire SED from ra¬ 
dio to y-rays allows one t o perform a reasonable estimation of 
these physical parameters (iTavecchio et al.l[T998l) . For the case 
of synchrotron emission alone, R, B and 6 will only decide the 
spectral normalisation. On the other hand, the shape of the ob¬ 
served spectrum is determined by the corresponding form of the 
underlying particle spectrum. To obtain further insight into the 
emitting particle distribution, we developed synchrotron emis¬ 
sion models with different particle distribution and incorporated 
them into XSPEC spectral fitting software. Particularly for this 
study, we consider the following particle distributions: 


(i) Simple power law (SPL): In this case, we assume the elec¬ 
tron distribution to be a simple power law with a sharp high 
energy cutoff, given by 


N{y)dy = Ky ^dy, y < y^ax (6) 

Here ymc^ is the energy of the emitting electron, p is the 
particle spectral index, K is the normalization and ymaxfnc^ 
the cut-off energy. Among these parameters, p and ymax are 
chosen as the free parameters. 

(ii) Cutoff power law (CPL): Here the underlying particle dis¬ 
tribution is assumed to be a power law with index p and an 
exponential cut-off above energy ymfnc^ given by 

N{y)dy - Ky^^exp(—)dy (7) 

ym 


For this distribution, p and y^ are chosen as the free pa¬ 
rameters. 

(iii) Broken power law (BPL): The particle distribution in this 
case is described by a broken power law with indices p\ 
and p 2 with a break at energy yi, given by 


N(y)dy — 


Ky P'dy, ymm<y<yb 

Kyf-P^^y-P^-dy, yh<y< y,„ax 


( 8 ) 


Here, p\, pi and yt are chosen as the free parameters. 

(iv) Log parabola (LP): For this case, the particle distribution 
is chosen to be a log parabola, given by 


N{y)dy = K{y (9) 


with Op and Pp chosen as the free parameters. 

We fitted the observed combined X-ray spectrum from Swift- 
XRT and MrSTAR with the synchrotron emission due to these 
different particle distributions as given above. A poor fit statis¬ 
tic with large reduced is encountered for the case of SPL 
since it fails to reproduce the smooth curvature seen in almost 
all spectral states (/I - /13). For CPL, the statistics improved 
for many states (Table O with lowest reduced of 1.01 during 
state D; whereas, the largest reduced is 1.41 for the state f6. 
The ht statistics improved considerably for many states for the 
case of BPL except for f6 and f8 corresponding to peak X-ray 
flare (TableO. In Figure|5j we show the cross plot distribution of 
the power law indices, p\ and p 2 , of BPL model during different 
spectral states. The index pi is poorly constrained for the state 
f3 due to the absence of Swift-XRT observation during this pe¬ 
riod; whereas, for state fll neither pi nor p 2 are well constrained 


due to the absence of AmSTAR observation. However, for most 
of the states, a strong correlation between pi and p 2 is observed 
(rs = 0.79 and p = 0.0018). Among all these particle distribu¬ 
tions, the best statistics is obtained for the case of the LP model 
with the reduced x^ decreased considerably during the flaring 
states, f5, f6 and f8 (TableO. Further, the reduction of one free 
parameter in case of LP with respect to BPL enforces the latter 
to be the most preferred particle distribution. In TablelH we give 
the best fit parameters for the case of the LP particle distribution. 
We see similar correlations as discussed in Section HTTl In fact, 
there is a strong linear correlation (p < le - 16) between the cor¬ 
responding parameters of the photon and the particle spectrum. 


5. Discussions and Conclusions 


In the present work, we performed a detailed study of the bright 
X-ray flare of Mkn 421 observed during April 2013 along with 
information available at other wavebands. We noticed the X-ray 
flare is not signihcantly correlated with the UV and, in addition, 
the variability amplitude of the former is considerably larger 
than the latter. This suggests that probably the X-ray and UV 
emission may belong to emission from different particle distri¬ 
butions. 

A detailed spectral analysis of the X-ray observations over 
different time periods during the flare suggests the emission to 
arise due to synchrotron mechanism from log parabola particle 
distribution. Though this particle distribution can be statistically 
more appealing, a broken power law particle spectrum cannot be 
excluded. Further, we extended the best fit LP and BPL particle 
distributions to low energies and predicted the UV synchrotron 
flux. During low X-ray flux states, the predicted UV flux agrees 
reasonably well with the observed flux. However, at high X-ray 
flux states, the observed UV flux is significantly higher that the 
predicted one by a factor of 2 - 50, with the larger deviations 
corresponding to the LP model. In addition, the variability of the 
predicted UV flux is much higher {Fvar,pred - 1.10 + 0.10) than 
the one obtained from the observed UV flux. This study again 
questions the similar origin of X-ray and UV emission. 

A plausible interpretation of this inconsistency between X- 
ray and UV fluxes can be done by associating the UV emission 
from the putative accretion disk. However, such thermal emis¬ 
sion from the disk w as never a importa nt contribution in the UV 
bands for Mkn 421 (lAbdo et alJl201 ill , and the UV spectral de¬ 
tail is not sufficient enough to assert this interpretation (Figure 
|6|l. Alternatively, the underlying particle distribution can be more 
complex than the ones studied in this work. Nevertheless, such 
a particle distribution demands a concave spectrum which is not 

J iQssible with our pre sent understanding of particle acceleration 
Sahavanathanll20()8l) . Hence, we attribute this unusual X-ray - 
UV behaviour of the source as a result of two popul ation electron 
distrib ution. A similar conclusion was obtained bv lAleksic et al.l 
(1201 4t) by studying a flare of the same source in March 2010. 
Such different electron distributions can be obtained if multiple 
emission regions are involved the emission process. If the flaring 
region is located at the recollimation zone of the jet, then a com¬ 
pact emission region ca n be achieved where the recollimation 
shock meets the jet axis (ITavecchio et alJl201 UlKushwaha et^ 


20141). Alternatively, episodic particle acceleration suggested by 


Perlman e t alJ (l20 06t) can be a reason for the second particle dis- 
tribution. fPerlman et al.l (l2005h have shown that this can also ex¬ 
plain the relative less variability of the optical/UV bands as com- 
par ed to the X-ray. 

iPian et aP (1201 4l) studied the same flare, starting from MJD 
56397, with emphasis on INTEGRAL and Termf-LAT data. They 
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also found trends of spectral hardening with flux. However, un¬ 
like our results, the INTEGRAL spectral data could be well fit by 
a broken power law spectral form. They modelled the broadband 
SED with a simple one zone SSC model, and required large vari¬ 
ations in the magnetic fields and Doppler factors to fit the SED 
of different states successfully. A lso, our results do not match 
with that of iMassaro et al.l (l2004l) as we do not see any correla¬ 
tion between as and /3j. Thus, our results are inconsistent with 
statistical particle acceleration. 

Our results indicate the potential of using broadband X-ray 
data to constrain the underlying particle spectrum and differenti¬ 
ate separate variable spectral components, especially if there are 
simultaneous data in other wavebands. Eigure [T] shows the en¬ 
ergy spectrum (in vFy) and the residuals for the BPL, the CPE 
and the LP for the state f6. While the residuals show structure 
as a function of energy for energies < 0.7keV, especially for 
the CPL and the BPL, the data and model agree within 10% 
at these energies. Moreover, we also implemented a fit at ener¬ 
gies < 0.7keV and verified that our results do not change sig¬ 
nificantly. Thus, this effect is unlikely to be a serious source of 
error in this work. More importantly, there is a clear systematic 
deviation from the data for the BPL and the CPL at the higher 
energies, beyond 5QkeV. The forthcoming satellite ASTROSAT 
(ISingh et al.ll2Q14li . with its wide band X-ray coverage and si¬ 
multaneous optical-UV measurement can be expected to make 
significant breakthroughs in this field. 
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Table 1: Details of Nustar pointings 


obsid 

start date and time exposure (in sec) 

60002023023 

2013-04-10 20:53:07 

118 

60002023024 

2013-04-10 21:26:07 

5758 

60002023025 

2013-04-11 01:01:07 

57509 

60002023026 

2013-04-12 20:11:07 

441 

60002023027 

2013-04-12 20:36:07 

7630 

60002023029 

2013-04-13 21:36:07 

16510 

60002023031 

2013-04-14 21:41:07 

15606 

60002023033 

2013-04-15 22:01:07 

17278 

60002023035 

2013-04-16 22:21:07 

20279 

60002023037 

2013-04-18 00:16:07 

17795 

60002023039 

2013-04-19 00:31:07 

15958 


Table 2; Details of Swift pointings. The exposure times are rounded off to the nearest seconds 


obsid 

start date and time 

xrt exposure 

uvot exposure 

bat exposure 

00035014061 

2013-04-10 02:04:58 

1079 

1066 

1085 

00080050016 

2013-04-11 00:30:59 

1118 

1076 

1128 

00032792001 

2013-04-11 03:41:30 

3488 

3468 

3502 

00080050017 

2013-04-11 21:48:59 

1449 

1419 

1453 

00080050018 

2013-04-12 00:33:59 

8726 

8635 

8746 

00080050019 

2013-04-12 21:53:58 

9546 

9428 

9572 

00032792002 

2013-04-14 00:38:59 

6327 

6253 

6362 

00035014063 

2013-04-14 23:04:59 

4942 

4874 

4965 

00035014062 

2013-04-15 23:07:59 

534 

522 

540 

00035014064 

2013-04-16 00:43:59 

10262 

10108 

10302 

00035014065 

2013-04-17 00:46:59 

8842 

8731 

8857 

00035014066 

2013-04-18 00:49:59 

6887 

6798 

6907 

00035014067 

2013-04-19 00:52:59 

6132 

6060 

6152 

00035014068 

2013-04-20 00:55:59 

5543 

5482 

4640 

00035014069 

2013-04-21 07:33:59 

394 

389 

397 


Table 3; Reduced values for the various time bins (fl to fl3), as marked in Figure 3, for the different models of the photon 
spectrum (broken power-law (bknpo) and log parabola (logpar)) and the particle spectrum (cutoff power law, broken power-law 
(bknpo) and logparabola) respectively. Observation ids of XRT and Nustar are also given. 


State 

1 

XRT 

Obs id 

Nustar 

Photon spectrum 
bknpo logpar 

Particle spectrum 
CPL BPL LP 

fl 

00080050016 

60002023024 

1.27 

1.02 

1.07 

1.06 

1.03 

f2 

00032792001 

60002023025 (a) 

1.40 

1.16 

1.26 

1.19 

1.16 

f3 

- 

60002023025 (b) 

0.97 

0.92 

1.01 

1.06 

1.01 

f4 

00080050017 

60002023025 (c) 

1.22 

1.04 

1.08 

1.07 

1.05 

f5 

00080050018 

60002023025 (d) 

2.33 

1.05 

1.21 

1.16 

1.06 

f6 

00080050019 

60002023027 

2.93 

1.08 

1.41 

1.39 

1.09 

f7 

00032792002 

60002023029 

2.12 

1.10 

1.21 

1.19 

1.10 

f8 

00035014063 

60002023031 

2.45 

1.14 

1.20 

1.31 

1.13 

f9 

00035014064 

60002023033 

1.78 

0.93 

1.06 

1.09 

0.93 

flO 

00035014065 

60002023035 

1.75 

1.10 

1.09 

1.10 

1.10 

fll 

00035014066 

60002023037 

1.23 

1.00 

1.03 

1.06 

1.00 

fl2 

00035014067 

60002023039 

1.14 

1.01 

1.08 

1.02 

1.01 

fl3 

00035014068 

- 

1.07 

1.07 

1.04 

1.07 

0.92 
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Table 4; Best fit parameters of the log parabolic photon and particle spectrum respectively. 


state 


Photon spectrum 

CYp 

Particle spectrum 
Pp 

P 

fl 

2.21 ± 0.02 

0.39 ± 0.01 

0.534 ± 0.027 

3.19 ±0.06 

1.96 ±0.11 

0.50 ± 0.04 

f2 

2.21 ± 0.01 

0.39 ± 0.01 

0.543 ± 0.023 

3.18 ±0.04 

1.98 ± 0.09 

0.50 ± 0.03 

f3 

2.14 ±0.05 

0.27 ± 0.03 

0.550 ± 0.074 

3.11 ±0.14 

1.25 ±0.16 

0.36 ± 0.06 

f4 

1.92 ±0.01 

0.38 ± 0.02 

1.254 ± 0.0630 

2.45 ± 0.06 

1.97 ±0.11 

0.77 ± 0.06 

f5 

1.82 ±0.01 

0.44 ± 0.01 

1.585 ± 0.0229 

2.10 ±0.02 

2.35 ± 0.04 

0.95 ± 0.03 

f6 

1.80 ±0.01 

0.46 ± 0.01 

1.634 ± 0.0205 

2.01 ± 0.02 

2.50 ± 0.03 

1.00 ±0.02 

f7 

2.08 ± 0.01 

0.41 ± 0.01 

0.805 ± 0.015 

2.82 ± 0.03 

2.10 ±0.04 

0.64 ± 0.02 

f8 

1.61 ±0.01 

0.37 ± 0.01 

3.272 ± 0.0668 

1.66 ±0.03 

1.92 ± 0.04 

1.23 ± 0.05 

f9 

1.99 ±0.01 

0.33 ± 0.01 

1.018 ± 0.0201 

2.69 ± 0.02 

1.62 ± 0.03 

0.61 ± 0.02 

flO 

1.83 ±0.01 

0.30 ± 0.01 

1.872 ±0.039 

2.33 ± 0.02 

1.46 ± 0.03 

0.77 ± 0.02 

fll 

2.33 ± 0.01 

0.30 ± 0.01 

0.280 ± 0.013 

3.55 ± 0.04 

1.44 ± 0.07 

0.29 ± 0.02 

fl2 

2.39 ± 0.01 

0.33 ± 0.01 

0.259 ± 0.012 

3.67 ± 0.04 

1.59 ± 0.08 

0.30 ± 0.02 

fl3 

2.25 ± 0.01 

0.35 ± 0.03 

0.445 ± 0.049 

3.30 ± 0.05 

1.78 ± 0.23 

0.43 ± 0.06 
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90 92 94 96 98 100 102 

MJD [56300+] 

Fig. 1: Multiwavelength lightcurve during MJD 56390 to 56403 showing : Panel 1: V band magnitude; Panel 2: UV flux in mJy; 
Panel 3: Degree and angle of optical polarization; Panel 4; Swift-XRT flux in counts/sec; Panel 5: MAXI flux in counts/sec; Panel 
6; AmSTAR flux in counts/sec; Panel 7; Swift-BA^ flux in counts/sec; Panel 8; Fenni-LAI flux in phjcnP'jsec. One Swift-XRT (or 
AmSTAR) point is plotted for each snapshot. 


90 

60 

30 

00 












Sinha et. al.: Study of Mkn 421 during huge X-ray flare in April 2013 



Fig. 2: Plot of hardness ratio as observed by AmSTAR FPMA. The upper panel shows the hardness ratio, and the lower one, the 
flux(counts/sec). There is a trend of spectral hardening with flux. 


1 2 3 4 5 6 7 8 9 10 11 12 13 



MJD [ 56300+] 

Fig. 3: The 13 time bins for which spectra have been extracted. 
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flux [ergs/cm /sec] 
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Fig. 4: Cross plot of parameters for the logparabolic photon (green stars) and particle (red circles) spectrum respectively. The peak 
curvature Ep clearly anti-correlates with alpha, and strongly correlates with the flux. However, there is no correlation between Ep 
and the curvature parameter /?, or between /? and a. 



Fig. 5; Cross plot of the 2 indices of the broken power law particle spectrum. The green points are not from combined spectrum, but 
are NuSTAR only and Swift-XKT only (see text). 
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Fig. 6: Observed X-ray and UV spectrum during two states. While the X-ray flux changes by more than a factor of 10, the UV flux 
remains constant. 



(a) cutoff power law, CPL 

Unfolded Spectrum 


(b) broken power law, BPL 



1 10 
Energy (keV) 


(c) log parabola, LP 


Fig. 7: Fitted spectrum and residuals for the three models for state f6. The differences start showing up above 50keV 
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